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Abstract. The behavior of annealed branched polymers near adsorbing surfaces 
plays a fundamental role in many biological and industrial processes. Most importantly 
single stranded RNA in solution tends to fold up and self-bind to form a highly 
branched structure. Using a mean field theory, we both perturbatively and numerically 
examine the adsorption of branched polymers on surfaces of several different geometries 
in a good solvent. Independent of the geometry of the wall, we observe that as 
branching density increases, surface tension decreases. However, we find a coupling 
between the branching density and curvature in that a further lowering of surface 
tension occurs when the wall curves towards the polymer, but the amount of lowering 
of surface tension decreases when the wall curves away from the polymer. We find 
that for branched polymers confined into spherical cavities, most of branch-points are 
located in the vicinity of the interior wall and the surface tension is minimized for a 
critical cavity radius. For branch polymers next to sinusoidal surfaces, we find that 
branch-points accumulate at the valleys while end-points on the peaks. 


1. Introduction 

Branched polymers play important roles in many biological and industrial systems, 
notable among them single stranded RNA (ssRNA) that in solution takes on a branched 
secondary structure [H121131 SI El EKHIHlEmo]. Recent experiments on viruses show that 
some viral RNAs, in particular, assume highly branched structures mm- The physics 
of polymer adsorption on different kinds of interfaces has, specifically, attracted a lot of 
interest for over half a century [131 [HI HSl HEl HZl EHl [HI EOl EH EH EH El] . In particular, 
it has been shown that polymer topology can effect the thermodynamic behavior of 
polymers near surfaces |25|- More recently the adsorption of RNA onto spherical gold 
nano-particles has been the focus of intense research because of its potential application 
in drug delivery or gene therapy [26l EH EHl EHl EHl EH EH EH El] • 

RNA is considered as an annea/ed branched polymer mainly due to the fact that the 
base-pair binding in RNA is often weak enough that the branching structure can change 
due to thermal fluctuations [HESj- Beyond the adsorption of RNA onto nano-particles, 
the behavior of annealed branched polymers next to surfaces of complex geometries 
is intriguing. Despite the presence of several excellent books and review articles, the 
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impact of branching on adsorption of biopolymers at planar or rough substrates is yet 
not well-studied. 

Several experiments compare the efficiency-directly connected to the free energy-of 
encapsidation of linear polymers and viral RNAs by virus capsid proteins [36l (STJ l38] . 
Field theoretic models have been used extensively to calculate the free energy of linear 
polymers [39l SOI HU HSl 031 011 05l 06l 07]. In a 1972 seminal paper de Gennes noted 
an equivalence between the statistics of a self-avoiding polymer and the n —)■ 0 limit of 
an 0{n) model of a magnet [IH], see Appendix A for a review of 0{n) model. Using 
this observation, which relates a mathematically interesting but unphysical limit for the 
model of a magnet to the statistics of polymers, it became possible to use the tools of 
statistical held theory to describe the physical properties of a polymer solution. Later, 
de Gennes held theory was expanded to describe the statistics of annealed branched 
polymers [IHl [50] . 

In this paper we use a mean held theory to study the adsorption of annealed 
branched polymers on diherent types of surfaces from a semi-dilute polymer solution in 
a good solvent. We study the ehect of curvature by examining the adsorption onto the 
exterior and interior of a spherical surface, and investigate the impact of roughness 
by examining the adsorption onto a sinusoidal surface. Instead of considering a 
random roughness, we employ a grating geometry because of its enormous mathematical 
simplihcations and the fact that it has been shown that qualitatively the essential 
features of the results are the same [23l 02] . 

By numerically solving the relevant nonlinear equations we hnd that compared 
to the adsorption to a hat wall, branching density, surface tension, and the monomer 
density all increase if the polymer is adsorbed onto the interior wall of a spherical cavity 
but decrease if adsorbed on the exterior surface of the sphere. While our results show 
that surface tension always decreases as branching density increases independent of the 
geometry of the wall, we hnd the interplay of curvature and branching density conspires 
to further lower the surface tension when the wall curves toward the polymers but 
lessens the amount of decrease in the surface tension when the wall curves away from the 
polymer. In the limit of large spheres, we solve the nonlinear equations perturbatively, 
which match the numerical results. Furthermore, in case of sinusoidal surfaces, we hnd 
inhomogeneity in the branching density as it increases in the valleys but decreases in 
the peaks. 

The remainder of this paper is organized as follows. In the hrst section we present 
our mean held approach and in the following section we use this method to investigate 
the impact of branching combined with surface curvature on polymer adsorption. In 
particular we will examine what ehect the branching structure has on the adsorption of 
polymers to nano-spheres and sinusoidal surfaces. We will hnish with a brief summary 
and present our main conclusions. 

In the appendix, for completeness and pedagogical reasons, we derive a simple 
held theoretical model for a branched polymer by revisiting the held theory developed 
by Isaacson and Lubenskv[49j for branched polymers and will spell out in detail the 
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2. Mean Field Approximation 

To describe a branched polymer on a lattice, we assume the polymer system consists of 
branch-points and end-points lying on the lattice sites, and bonds that join neighboring 
lattice sites. We treat the system as an annealed branched polymer, so the structure of 
the branched polymer is not hxed. For simplicity, we assume that all branch-points are 
exactly of order three because all higher order branch-points can be considered as many 
order three branch-points in close proximity to each other. For example, two order three 
branch-points close together will behave very similarly to an order four branch-point. 
The quantities that describe such a polymer system are: (i) Np, the number of polymers; 
(ii) Nfj, the number of bonds; (hi) Ni, the number of endpoints; (iv) N^, the number of 
branch-points; and (v) Ni, the number of loops. There is a constraint [19] relating most 
of these quantities such that 

2{Np - Ni) = W - Ns. (1) 

In this paper, we consider a system of branched polymers with no loops and set Ni = 0. 

The primary statistical quantity of interest is the multiplicity fl{Nb, Ni, Ns;V), 
dehned as the number of ways to arrange a polymer system of W bonds, Ni end-points, 
and Ns branch-points on a lattice that occupies a volume V. This quantity is equivalent 
to the number of microstates for the microcanonical ensemble. From the multiplicty, we 
can form the grand canonical partition function 

E(K,f„h-,V) = 

'^I<'^‘fPf[‘n(Ni„N,,N3,V), ( 2 ) 

Nt,Ni,N3 

where K, fi, and fs are the fngacities for the bonds, end-points, and branch-points 
respectively. Note that we use Eq. along with the assnmption that A) = 0 to 
eliminate the dependence npon Np. From this dehnition it is simple to derive expressions 
for the nnmber of bonds, end-points, and branch-points as derivatives of the logarithm 
of the grandpartition fnnction 

Ni, = KNinE{K,fuh), (3) 

iVi=/ij((rliiH(A'./i,/3). (4) 

Cjl 

iV3 = /3j^liiH(A'./i,/3). (5) 

ofs 

Following the idea of de’Gennes |1H], and using the methods of Lubensky and Isaacson 
[l9], we equate the grand partition function for the branched polymers system with the 
n —)■ 0 limit of the partition fnnction of an 0{n) model of a magnet 

E{KJ,Js;V)^Z{K,fsJs;V). 


( 6 ) 
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The partition for the 0{n) model of a magnet can be written as a function integral over 
a continuous field 

Z{K,hJ,;V) = jvf; (7) 

where 77 is the effective Hamiltonian 

X ' x' 

- ylV’(®)r + 

2 

- fi\fKfji{x) - . (8) 

We emphasize here that the parameters K, fi and /a take on different meanings for the 
0{n) model of the magnet (see Appendix A). For example, the K parameter represents 
the coupling constant between the nearest neighbors spins in the 0{n) model of a magnet 
and the ^|^ held is the average magnetization in a small region for the magnet. However, 
the t/? held in Eq. (j8]) is proportional to the monomer density for the branched polymer. 
The hrst term in Eq. (0) is an entropic term that smoothes out the held. The 

in the hrst term is the inverse of the nearest neighbor operator S(^x,x'), dehned as 
1 if x and x' are neighboring sites on a lattice and 0 otherwise. The second term in 
Eq. (IH]) proportional to |i/’P is due to the nearest neighbor attraction, which tends to 
increase the ^/^(x) held. The third term proportional to is repulsive representing 
the self-avoiding nature of the polymer, and leads to a decrease in the 'ij){x) held. The 
fourth and hfth terms proportional to 'ipi{x) and ipfix) are both attractive, and show 
that both branch-points and endpoints serve to increase the 'if{x) held. It is the balance 
of these attractive and repulsive terms that creates a well-dehned hnite 'ijj{x) held in 
equilibrium. For completeness, as well as pedagogical reasons, we show the derivation 
of the equivalence between the grand canonical partition function, Eq. [21 for the polymer 
system and the partition function for the n —)■ 0 limit of the 0{n) model and all relevant 
approximations for hnding the ehective Hamiltonian, Eq. [HI in Appendix A. 

We now make a mean held approximation and assume that the value of the held 
'ip{x) is uniform and is well approximated by its average value i/jq. Thus the sum over 
the nearest neighbors simply becomes 

V 

'^5{x;x')'iIj{x') = Z'ljjQ, (9) 

X 

with 2 ; the number of nearest neighbors or coordination number. The inverse of the 
nearest neighbor delta function is then simply the reciprocal of the coordination number 



Adsorption of annealed branched polymers on curved surfaces 


5 


El = So in the mean field theory, the grand canonical partition function 

is the exponential of the effective Hamiltonian evaluated at its minimum 

E{K, /i, /s; H) ^ (10) 

with -^0 found by minimizing Eq. (E]), = 0, 

(t^Ky,-h•/K-h^^Pl + !A,p| = o. ( 11 ) 

From stastical mechanics we can identify the grand potential, (3^g{K, fi, f^;]/) = 
-ln[E{K, fi, f^;]/)] using Eqs. ([H]), dSD and ([10]) , 


P<l>G{K,f,,fs;V) = 




- Kji>l - VVKiPo 


hKi ..K\ 
—V-o + 


. ( 12 ) 


with a the lattice spacing. Inserting Eq. flT^ in Eq. (|2|) and using Eq. (ITT]) we find the 
average monomer density 

Ci = N,/V = (13) 

At this point it is convenient to define a new field 4>{x) such that in the mean field 
approximation the average value 0 q = 'iI)q/2zo? is the monomer density. The grand 
potential written in terms of the new field is 


/3«I>g(J^,Fi,F3;E) _ 

V 

(1 - zK)(f^ - /i\/^^0o 

V CL 

- ^f3{2zaK)^(pl + ^K^z'^a^cpl- (14) 

By comparing Eq. (lT4j) with the expression for the grand potential for a linear polymer 
in good solvent |12], we can identify r = (1 — zK) as the chemical potential of 
monomers (such that djd^G/df = Cft) and v = K^z^a^ as the excluded volume. It 
is also convenient to absorb \/2zK and { 2 zK )2 constants in the end and branch point 
fugacities, respectively such that the grand potential can be written in a much simpler 
form 

(d^G{KJij3;V) 

V 

+ (15) 

v O'* t) z 
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The parameters /i and /s are to physical quantities. Using Eqs. dl]) and ([S]), we find 
the end-point and branch-point densities 

Cl = ^00, (16) 

Va 

C3 = ^/sV^^o- (17) 

3. Results and Discussion: Branched polymers adsorption onto different 
surfaces 

We now apply the field theory presented in the previous section to a semi-dilute system 
of annealed branched polymers and investigate their adsorption to different surfaces. 
More specifically, we consider a solution of branched polymers with a monomer density 
Cb, where the polymers all have a fixed length L, and a tunable average branching 
number W- 

The adsorption mean field energy of the branched polymer, F — Fq can then be 
written as 


F-Fo = --fa^ j 

The first term in Eq. flTHD is a surface integral that gives the contact energy between the 
surface and the polymer. The first term in the volume integral is associated with the 
entropic cost of a non-uniform polymer distribution [51]. The rest of the terms in Eq. 0181 ) 
are related to the free energy of a branched polymer in mean field approximation, see 
Eq. ([I5|). 

Considering the constraint that the total number of monomers is fixed 

j dV = constant = N, ( 19 ) 

Eq. ffTSj) can be rewritten as 
F-Fo = j dS(f^+ 

- ~ + /3^(0^ - 4^^)) - A(0^ - Cb)^ (20) 

where A is the Lagrange multiplier. 
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Minimizing Eq. fl20p with respect to the held 4>{x) gives the following Euler- 
Lagrange equation 

= -A0 + 1^0^ - -^(/i + \h4>^) (21) 

6 2\/af‘ 2 

subject to the boundary condition 


[h.V0 + (6007)0]^ = 0 


( 22 ) 


Eor simplicity, we rescale the Lagrange multiplier E = ^ and introduce a length 
that characterizes the strength of the attraction between the surface and monomers 
as K~^ = The other boundary condition is natural, far from the surface the held 

should be uniform, V^0 —?■ 0, and and take on the bulk value 0(x) yjcf. Using Eq. [2ll 
the Lagrange multiplier can be written as 


E 



2 Cfe 


2cJ 


(23) 


where ci = f\c^‘^l\fcA and C 3 = /6 are respectively the end-point and 

branch-point concentrations far from the surface. For long polymers with no loops 
in which Ni = + 2, the second and third terms in Eq. fl23l) correspond to the 

ratio of the end-points and branch-points to monomers. To make all the quantities 
dimensionless, we rescale the held with respect to the bulk value 9{x) = 0 (x)/yTb and 
the spatial coordinate with respect to the Edwards correlation length |S2] x = x^e, where 
^E = 7^^5 the equation of motion simply becomes 


= 2{e^ -6) + Ai{e -i) + sA^ie - e^) 


(24) 


with Ai = (ci/cfe)/(z/Cb), A 3 = {c 3 /cb)/{i^cb). 

The Ai and A 3 quantities measure the relative importance of the branching 
structure of the polymer to the steric ehect or excluded volume interaction between 
monomers in solution . The numerators ci/cb and cs/cf, are the ratios of the 
concentrations of end-points and branch-points to the total number of monomers, 
respectively. For large polymers, these ratios can approach 1/2 for a maximally branched 
polymer. The denominator vcb is a hlling fraction of the polymer; dilute solutions will 
have small values of vCb, and a dense polymer system with no solvent will have a value 
of one. With the new scaled coordinates, the boundary conditions become 

= 0 (25) 

lim 9 = 1 (26) 

tc—)-oo 

with K = k^e- In terms of the the new mean-held order parameter, 9, we can rewrite 


89 

dh 


k9 
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F - Fo = ^fECb(^-H j dSe^ 

+ j dv{(yef + {6^ - if 

- Ai{e - 1 ) - A3{2e^ - 39 + 1 ))^ (27) 

In the following sections we employ Eq. fl27|) to calculate the free energy of a polymer 
next to a flat, curved, spherical and sinusoidal surfaces. We note that while the Lagrange 
multiplier A acts like a chemical potential in open surfaces for hxing the density of bulk 
polymers, in the closed systems, like inside a sphere, it’s used to fix the number of 
monomers inside the shell. 

3 . 1 . Analytical Calculations 

3 . 1 . 1 . Flat wall Next to a flat wall, the Euler-Lagrange equation, Eq. fl2^ . subject 
to the boundary conditions given in Eqs. fl25|) and fl26|) can be solved perturbatively. 
We assume the attractive interaction between monomers and wall is smaller than the 
monomer-monomer repulsion (k 1). The solution to Eq. (124|) to the second order in 
K can then be written as 


0^1 + j{l + j)e-^^ (28) 

where A = \/4: + Ai — 3A^ with Ai and A^ proportional to the number of end (W) and 
branch {Nf) points as given below Eq. fl2T)) . For a single long polymer with no loops 
Eq. yields Ni = N 3 + 2. So we can simply write Ai - 3 A 3 = 21/(1 - Nf/iuN^) 
implying that Ai — 3A 3 < 0 and \Ai — 3 A 3 I < 4 should be satished for real solutions. 
For Ai = A 3 = 0, Eq. fl28l) converges to the profile of a linear polymer next to the flat 
wall [42]. As clearly shown in Eq. (j28|) . the density of branched polymers are larger than 
the linear ones next to a flat wall. 

Inserting Eq. fl28|) into Eq. (|27|) . we find the change in surface tension, energy per 
unit area, to the second order in k 

(29) 

where T^ is the difference in tension from a linear chain due to branching and is given 
by 


Tb{n, Ai, A3) 


Ai -|- 3 A 3 

A2 


)k 


,3A2-A3-4 + 2Ai + 12A3. 
' 2A3 




( 30 ) 
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Since the quantity Ff, for all acceptable values of Ai and A 3 is always positive, according 
to Eq. fl29l) the surface tension due to adsorption of a linear polymer to a flat wall is 
always higher than that of a branched one. In the next section, we calculate the impact 
of curvature on the polymer density profile and the free energy of the system. 


3.1.2. Curved wall To investigate the effect of curvature analytically, we assume that 
the flat wall is slightly bent to form a large sphere. The curvature could be either toward 
or away from the polymer. The radius of curvature b is considered to be large compared 
to the correlation length (6 = b/^E S> 1). We can obtain the perturbative solutions 
in spherical coordinates through Eq. [2l] assuming that 0 = 1 + 5 (5 <C 1) at the weak 
adsorption limit (k 1). According to the direction of wall curvature, polymer solution 
is considered either to be in the interior (in) or the exterior (out) of a sphere of radius 
b. The perturbative solutions are then 

kb 


Oout{r) ~ 1 + 


1 + Ab — kb 


(f) 


o-A{r-b) 


(31) 


for {h < r < 00 ) and 


Oi. 


i + sKi)- 


sinh Ar 


-^ (32) 

Abcosh{Ab) — sinh(Afe)(l + kb) 

for (0 < r < 6). Asumming 6 S> 1, we can write the monomer concentration on the 
surface to the second order in k 


e, 


out 


Q. ~ 


K 

A’ 

k 


ft, . 1 ft,s 

1 H—T (1-^ H—t) 

4V aJ 


K ■ 

A' 
k • 


Ab A’ 


(33) 

(34) 


Comparison of Eqs. fl33ll and fIMll with Eq. [2H] reveals that branched polymer 
concentration in the vicinity of a flat wall increases if the wall bends toward the polymer 
and decreases if the wall bends away from the polymer, consistent with the results 
obtained for linear polymers |12]. In order to obtain the change in surface tension due 
to the wall curvature, we insert the concentration profiles given in Eqs. fldTll and 
into Eq. fl2711 and keep terms up to the hrst order in 1/b. Then we have 


(d bo) out 






and 


(7 - 7o)i 


a^Cf, ^ / r . 


( 36 ) 
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Concentration 

Surface Tension 

Flat surface 

^'flat = 1 + + i) 

- S - f - y) 

Convex surface 

^out — ^flat ~ 

Ayout ~ Ayflat + + Tg) 

Concave snrface 

(^in = ^flat + ^ 

^7in ~ Ayflat g^^^ ( + Tg) 


Table 1. Summary of the analytic results. The table shows the impact of the wall 
geometry and the topology of the polymer (embedded in Fg and Ff,, see the text) for 
two physical quantities of 6 and Ay. The quantity 9 is related to the concentration of 
the polymer at the surface by ct,d^ with Cb monomer density. Ay is the change in the 
surface tension due to the presence of the wall. 


where 7 — 70 is the surface tension for the polymer next to a flat wall based on Eq. fl2^ 

- 2 

and ^ is the difference in surface tension due solely to the geometry of the wall. The 
quantity 


Tg{K, b, Ai, A^) 


8^2 -A^-16 + 4Ai + 36^3 

4M ^ 


(37) 


is due to the the coupling between the geometry and branched structure of the polymers. 
If we set Ai = 0 and ^3 = 0 then A = 2 and Eqs. fl3^ and fl36l) reveal the impact of 
wall curvature on the surface tension for the case of linear polymers. As expected, the 
surface tension decreases if the wall bends toward the polymer and it increases if the 
wall bends away from the polymer. 

Quite interestingly, the Eg expression shows the importance of coupling between 
wall curvature and polymer branching on the surface tension. A glance through Table 
1 . shows that the sum or difference Tf, ± Tg corresponds to the change in the surface 
tension due to the coupling between branching and wall curvature. Since this term 
is positive for both the convex and concave surfaces, we see that branching always 
decreases the surface tension. However, the coupling between branching and curvature 
further lowers the surface tension when the wall curves toward the polymers compared 
to when the wall curves away from the polymer. 


3.2. Numerical Calculations 

The above analytical calculations were related to the surfaces with large radius of 
curvature. To study polymer adsorption on the surfaces with higher curvature, we 
need to numerically solve Eq. fl2T|) for both the interior and exterior of smaller 
spheres. Since this is mainly a phenomenological model, we follow similar parameters 
to previous numerical works and values that are typical for virus capsids and RNA sizes 

[laEHlnasal ESIES EH]. 

3.2.1. Outside the sphere In this section, we consider smaller spheres and obtain the 
polymer concentration prohle outside the sphere vs the scaled distance from the surface 
of the sphere, r — b. We obtain numerical results by solving the nonlinear differential 
















Adsorption of annealed branched polymers on curved surfaces 


11 



Figure 1. Profiles of the scaled polymer density amplitude 6 vs scaled distance 
from the surface of the sphere r — b for several relative branching density A^. The 
solid line gives the profile for a linear polymer with no branching, and dashed lines 
give the profiles for branched polymers with branching density = 0.2,0.4,0.8,1.0. 
Other parameters used are k = 1, 6 = 5.0 and v = 0.5. The inset of Fig. [T] is the 
concentration profile for several values of the radius of the sphere (6 = 0.1,1,5,100) 
for the branching density A^ = 0.8 with k = 1 and v = 0.5. As the radius of the sphere 
increases, monomer concentration on the surface increases. 


equation Eq. fl2^ subject to the boundary conditions given in Eqs. fl25|) and fl26|) . The 
results are presented in Fig. [H which illustrates that both the monomer concentration 
at the wall and the thickness of the adsorption layer increase as the branching density 
increases. The surface excess adsorbed onto the sphere can also be calculated using the 
concentration prohle 9, 


fOO ~ 

nex = ^ECb dr{e‘^-l){jf (38) 

Jb 0 

Figure [ 2 ] shows the surface excess as a function of the branching density. As 
the branching density increases, the surface excess also increases. For = 

0.2, 0.4, 0.6, 0.8,1.0, the ratio of surface excess of a branched to a linear polymer is 
^branched!^linear =x.ll 8 , 1.264, 1.451, 1.696, 2.030, respectively. Inset of Fig. [2] is a plot 
of surface excess vs the radius of the sphere for different branching densities A 3 = 0 . 2 , 
0.4, 0.6, 0.8, 1.0. According to the inset of Fig. O for a given branching density, there is a 
critical radius r* for which the surface excess has a maximum. We hnd that the position 
of the critical radius decreases linearly as branching density goes up, i.e., r* oc —A 3 . 

Using the numerical solution for 6 and Eq. 0271) . the surface tension can be written 
as 

poo 

+ / dr((02-l)2_yl^(0_i) 

Jb 

- A3{29^ - 30 + 1) + (V0)2) . (39) 
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^^3 


Figure 2. Surface excess (in units of ^ECb) vs branching density for k = 1.0, b = 5.0 
and u = 0.5. Inset shows the excess on different size of spheres for k = 1.0 and u = 0.5. 
The solid line represents the linear chain with = 0, and the dashed lines give the 
surface excess for branched polymers with branching density = 0.2,0.4,0.6,0.8,1.0 
as dashing gets smaller. 



^3 

Figure 3. Surface tension (in units of ) vs branching density for k = 1.0, b = 5.0 
and V = 0.5. Inset is a plot of surface tension as a function of radius of the sphere 
for it = 1.0 and v = 0.5. The solid line represents the linear chain with ^3 = 0 , and 
the dashed lines give the surface tension for branched polymers with branching density 
A 3 = 0.2, 0.4, 0.6,0.8,1.0 as dashing gets smaller. 


Figure [3] is a plot of the surface tension as a function of the branching density, which 
shows that as the branching density increases, the surface tension decreases. The inset 
shows the tension vs the sphere radius for different branching densities A 3 =0.2, 0.4, 
0.6, 0.8, 1.0. As the radius of sphere increases, the tension decreases consistent with the 
perturbative results. 

3.2.2. Inside a sphere In this section, we obtain the polymer concentration prohle 
inside a sphere (r < b) by solving the nonlinear differential equation Eq. (]2T|) subject 
to the boundary conditions given in Eq. (12^ . In addition, because the polymer is now 
conhned inside an impermeable sphere, the total number of monomers, N is hxed. In 
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Figure 4. Concentration profile = 9(x)y/cb) inside the sphere for branching 

density = 0.05 (solid line) and = 1.04 (dashed lines) for 6 = 12, it = 1.0, z/ = 0.5 
and N = 3000 where the lengths are in units of a, the Kuhn length. Inset shows the 
details of the plot next to the surface. 



Figure 5. Surface tension for polymer adsorbed in the interior of the sphere at 
As = 0.05 (solid line) and A 3 = 1.04 (dashed lines) for R = 1.0, v = 0.5 and N = 3000 
where the lengths are in units of a, the Kuhn length and the energy is in units of fcsT. 


terms of normalized length scale and the order parameter 6 , we have 


N = / Anr'^O'^dr. 

Jo 

The snrface tension in this case can be written as 
7m - 7o = 77 ^ ( - 


+ / dr{{e^-l)^-A3{e-l) 


>0 


713(20^- 30 + 1 ) + (W)^)(-)^ 

b 


( 40 ) 


( 41 ) 


The concentration prohle as a function of r, the distance from the center of the 
sphere is shown in Fig. 0] for a branched polymer with branching density A 3 = 0.05 
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(solid lines) and A 3 = 1.04 (dashed lines). As illustrated in the hgure, due to the 
attraction between the wall and the polymer, more monomers are attracted to the 
surface of sphere as the branching density increases. 

Figure |5] illustrates the surface tension as a function of the radius of the sphere. 
Quite interestingly, we hnd that there is an optimal size for the radius of sphere for 
a given hxed chain length. For the parameters k = 1.0, u = 0.5 and N = 3000, the 
optimal radius is b = 11.25 with A 3 = 0.05 and is 6 = 10.90 for A 3 = 1.04. When the 
polymer is more branched, the optimal radius becomes smaller. This is mainly due to 
lower cost for the excluded volume interaction. 

3.2.3. Sinusoidal grating In this section, we consider a polymer solution next to a 
sinusoidal surface, 2 = zqcos{{2'k/ X)x). Here zq is the amplitude and A is the wavelength 
of the surface. As mentioned in the introduction, this should give some insight into the 
behavior of branched polymer next to a rough random surface because the qualitative 
features of the results are the same. To obtain the concentration prohle, 9 we solve 
Eq. fl24l) subject to the boundary condition given in Eq. fl25ll using a hnite element 
method in 2 D. The numerical results are shown in Figs. El [ 7 ] and [HI as contour plots of 
the polymer density next to the sinusoidal adsorbing surface. In all cases we keep the 
strength of the attractive interaction between the surface and monomers the same. 

Figure El shows that the prohle of a branched polymer next to a sinusoidal surface 
is similar to that of a hat wall if the wavelength of surface huctuations is large, for 
example X = X/^e = 20. The hgure also illustrates that the concentration of genome 
is higher in the valley compared to the peaks. This is consistent with the perturbative 
results presented in previous section, in that if the wall curves away from the genome, 
the monomer concentration decreases, otherwise, it increases. The non-uniformity in 
the concentration prohle at the wall becomes more apparent as we decrease A; i.e., the 
genome concentration becomes much higher at the valley compared to the peak, see 
Fig. [71 In the hgure, the amplitude of surface huctuations, zq = z/^e = 0.5, is chosen 
to be relatively small to emphasize on the diherence between the genome prohle next 
to the hat and sinusoidal walls. Note that the amplitude Zq = 0.5 in Fig. [ 7 ] is 10 times 
smaller than that in Fig. El Nevertheless, since A = 2 is 10 times smaller in Fig. d the 
impact of surface huctuations are more pronounced. 

In addition, we hnd that not only the concentration prohle at the wall is not uniform, 
the distribution of branch-points is not homogeneous either. Figure El illustrates the ratio 
of branch density to the monomer density at a sinusoidal surface. The hgure shows that 
the branching concentration is higher at the valley with respect to peaks consistent 
with our perturbative results in section 2 , where we found that the branching density 
increases if the surface is curved toward the polymer and decreases if the surface is 
curved away. 
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Figure 6. Contour plot of the concentration profile for the polymer with branching 
density A 3 = 0.8 for it = 1, A = 20 and io = 5. The solid thick black line is the position 
of the surface. The upper contour corresponds to a local concentration of 0^ = 1.25 
that increases to 6 ^ = 3.6 towards the bottom of the valley. The x and z coordinates 
are distances along and perpendicular to the corrugations respectively. 



X 

Figure 7. Contour plot of the concentration profile for the polymer with branching 
density A 3 = 0.8 for k = 1, A = 2 and zq = 0.5. The solid thick black line is the position 
of the surface. The upper curve corresponds to a local concentration of 0^ = 1.5 that 
increases to 0^ = 5 towards the bottom of the valley. The x and 5 coordinates are 
distances along and perpendicular to the corrugations respectively. 


4. Conclusion 

In this paper, we use the held theory methods based on the n —)■ 0 limit of an 0{n) 
model to describe the statistics of an annealed branched polymer. We, in particular, 
carefully examine the behavior of branched polymers next to various adsorbing walls 
both analytically and numerically. 

We show that for the annealed branched polymers, increasing the branching density 
will increase the concentration of polymer but decrease the surface tension next to hat, 
inward curving, and outward curving walls. In comparison to a hat adsorbing wall, we 
hnd that when the wall curves toward the polymer solution, the tension decreases but 
when it curves away from it, the tension increases. While these results are consistent 
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Figure 8. Contour plot of the ratio of branch to monomer concentrations with 
branching density = 0.8 for k = 1, A = 2 and zq = 0.5. The solid thick black 
line is the position of the surface. The upper curve corresponds relative concentration 
of branch points of ca/c = 0.4 and it increases to ca/c = 0.9 towards the bottom 
of the valley. The x and z coordinates are distances along and perpendicular to the 
corrugations respectively. 


with those found for linear polymers next to different type of walls, we interestingly 
found a correlation between the branching and curvature that causes a further lowering 
of surface tension when the wall curves towards the polymer, but decreases the amount 
of lowering of surface tension when the wall curves away from the polymer. 

Our numerical solutions for the adsorption of polymer to the exterior of small 
spheres show that increasing the branching lead to an increase in the surface excess. 
For the polymers adsorbed in the interior of small spheres, we find a minimum in the 
surface tension as a function of radius of sphere. This clearly demonstrates the interplay 
between monomer-wall attraction and monomer-monomer excluded volume interaction. 
Our hndings also indicate that branching decreases the optimal radius of sphere, as more 
polymers can sit in the vicinity of the wall without a huge cost for monomer-monomer 
repulsion. This result has a considerable consequence for the encapsulation of RNA 
by virus shell proteins, and could suggest a non-specific mechanism for the preferential 
packaging of viral RNA to cellular RNA in vivo [56]. 

Furthermore, we hnd similar effect for the adsorption of polymers onto sinusoidal 
surfaces. The concentration of polymers increases in the valleys compared to peaks and 
also branching density goes up in the valley section compared to the peaks. This effect 
becomes more pronounced as wave-length decreases. 

Understanding the mechanisms involved with the adsorption of annealed branched 
polymers onto different surfaces will play a critical role in biomedical technologies. In 
particular, the paper was inspired by the idea of using functionalized gold nano-particles 
to bind RNA for gene delivery [27], which has industrial applications for biosensors 
and microfluidic devices, and even possible medical application for gold nano-particles 
encapsulated by virus coats as potential tools for gene therapy. [Ml ETj IMl EHl IMl 1 ^ . 
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Appendix A: 0{n) model of a magnet 

In this appendix, we derive the equivalence between the grand canonical partition 
function, Eq. O for the polymer system and the partition function for the n —)■ 0 limit 
of an 0{n) model. Note that the 0{n) model corresponds to a magnet whose magnetic 
dipole of its atoms has n components . The Ising model commonly studied in most 
statical mechanics courses is the 0(n) model for n=l. The n=0 limit is an interesting 
case as it reproduces the statistics of a self-avoiding linear polymer. 

The Hamiltonian of the 0{n) model of a magnet is 

V 

H({5}, A, /i, h- V) = -KJ2 s.- S,, 

{x,x') 

V V 

(A,l) 

X X 

where Sx is an n dimensional vector of hxed length at each lattice point x and A, /i, 
and /a are the coupling constants. The hrst sum in Eq. lA.ll is over all pairs of nearest 
neighbors and and J 3 [S'a;] are the source terms for end-points and branch-points 

respectively. We have used some prescience in giving the coupling constants A, /i, and 
/s the same symbol as the fugacities in Eq. (Ej). The partition function for the 0{n) 
model is then 

^n(A,/i,/3;E) = (A.2) 

X 

where trg dehned as 


n(/ 

fs e"'® = ^- (A.3) 

n( / Sf - n)) 

i —oo i 


and the size of the spin is subject to the condition l^l = ^/n or Yh Sf 

_ i 

Using the power series dehnition of the exponential, Eq. (1A.2I1 can be 


= n. 

written as 


Z„(A,/i,/3;U) = 


OO 


E 

Nf,=0 


'W 


oo 


E 

N3=0 


xNs 

Is 


A 3 ! 


OO 


E 

Ni=0 


fNi 


Ai;U), 


(A.4) 
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Figure Al. Two possible configurations of 4 bonds, 3 end-points, and 1 branch-point 
on a 3 X 3 two-dimensional square lattice. The left diagram is not a valid conhguration 
while the right diagram is a valid configuration and is counted. 


with In defined as 


In{m,N^,Np,V) = 


V 


X 

V 

J2j3[s. 


V ^ N, 

5^ s. ■ s 

{x,x') 

Ns/ V ^Ni 

E MS 


V 


ntrs.W]V3!V,!Ecq{S}], (A.5) 


The comparison of the grand canonical partition functions defined in Eq. ([2]) with the 
partition function for the 0(n) model in Eq. flA.dp reveals the similarity between the 
two models. It is now obvious why the coupling constants in the 0(n) model were 
chosen to be labeled the same as the fugacities in the polymer system on a lattice with 
the excluded volume interaction. To show full equivalence between the grand canonical 
ensemble of self-avoiding branched polymers 5 in Eq. (E]) and the n —)■ 0 limit of partition 
function for the 0(n) model of a magnet (Eq. flA.2p . we only need to show that in the 
n —?■ 0 limit the expression gives the multiplicity hi, or counts the number of ways of 
arranging a self-avoiding branched polymer on the lattice. 

Each Cq,[{S'}] in Eq. flA.Sj) can be represented graphically. For the lattice, a product 
of neighboring S-vectors that lie on points x and x' can be represented by a line drawn 
between the points. The and source terms at point x and x' can be 

represented by circles and triangles placed on their respective points. 

As an example we present two possible configurations, one that represents a valid 
branched polymer configuration and one that does not, with Nb = 4, Ni = 3, and 
TVs = 1 on a 3 X 3 square lattice in Fig. lAll The corresponding Cq,[{/S'}] terms for both 
graphs in Fig. lAl I contain 4 dot products of neighboring S'-vectors (one for each bond), 
3 Ji[S] terms (3 circles), and 1 JslS] term (1 triangle). The Cq,[{S}] term for the right 
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C,1S] = (S, . s,) (S 4 ■ s,) (S, . Ss) 


(Ss ■ S,) J1IS2] Ji[S 4 ] JilSJJsISs]. (A.6) 


where the indices indicate which lattice site is associated with each term. 

We will show below that in the n —)■ 0 limit the trace of a single conhguration gives 
one if the graph corresponds to a physically valid branched polymer conhguration and 
zero otherwise, i.e., 


V 




1 for a valid, 

0 for a invalid. 


(A.7) 


This indicates, in the n —)■ 0 limit, the term given in Eq. flA.4l) counts the number 
of valid physical conhgurations, which is how multiplicity is dehned in Eq. ([2]). The 
condition in Eq. flA.7|) establishes the equivalence between the partition function for 
the 0(n) model in Eq. (lA.4p with the grand canonical partition function for a hexible 
branched polymer in Eq. 1 


t.e. 


A, A; E) = S(K, A, A; E). 


(A.8) 


Evaluation of the trace 

The trace of the conhgurations Cq,[{S'}] in Eq. flA.7j) takes the form of products of S- 
vectors on each lattice site x. These products can be evaluated by using the following 
generating function 


trs Fii • • • S'. = 


The trace of the generating function tr^ can be evaluated in closed form 


d 

dki^ 

kS 


d 


dki 


trs e’^'^ 


(A.9) 


fc ^.0 


CO 

tr r'^S _Sr' ^2. 


1 / n\k\ 


^ r(i 

p=0 ' ■ 


p) p\ 


2\P 


(A.IO) 


In the limit of n —?■ 0 it simplihes to the form 

limtrse'"^ = 1 +-IfcP. (A.ll) 

n^O 2 

A detailed step by step derivation of Eq. flA.lOD is presented in appendix A of Ref. [57] . 
Inserting Eq. (lA.lR into Eq. (1A.9I1 . we hnd the product of the components of the S- 
vectors in the absence of source terms is 


lim trs S'ii ■ • • S'. 

n^O 


1 if p = 0, 

if P = 2, 

0 otherwise. 


(A.12) 
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The expression in Eq. (IA.12P evaluates to a non-zero value only in the lattice sites with 
exactly 0 or 2 bonds terminating on them. 

To describe the generating function for end- and branch-points, we construct Ji 
and J3 functions, respectively, such that they satisfy the following equations 

lim tr^ S*!Tils'] = lim Sij, 


limtrsSiSjSfcJ3[S] = lim y/riy^ Sudjidku 

n^O n—s-fl .i—' 


n^O 


(A.13a) 

(A.13b) 


while all other traces involving the sources such as tr^ SjSj Ji[S], tr5SjJ3[S] and 
trg SjSj J3[S] are equal to zero in the n —)■ 0 limit. Using Eq. flA.lOp it is straightforward 
to derive the following expressions for Ji and J3, 


i 

A ^ 


AS = 


3n2 


-|S'2 

n 


S + 3S). 


(A.14a) 


(A.14b) 


The sum in Eq. flA.13aP evaluates to a non-zero value only if for every lattice site 
containing the source term Ji there exists exactly one bond terminating on that site. 
Similarly the sum in Eq. flA.13bp evaluates to a non-zero value only if for the sites with 
branch-point J3, there are exactly three bonds terminating on that site. The factors of 
1/y^ in Eq. flA.13aP and ^/n in Eq. flA.13bp enforce the no-loops condition necessary 
to obtain the statistics of a self-avoiding polymer. 

In general, for every valid conhguration with Np connected graphs, A^i end-points, 
and N 3 branch-points the product of all the lattice site integrals gives 


V 


limTTUsC„[{5}l = 


n—>0 


(A.15) 


Using Eq. (ID. the exponent in Eq. flA.lSp is equal to the number of loops, so in the 
n —)■ 0 limit only those valid conhgurations with no loops evaluate to 1, while all others 
evaluate to 0. 

It is important to note that the conditions given in Eqs. flA.13aP and flA.13bp do 
not forbid the multiple source terms sharing the same lattice site. This oversight can 
be remedied by redehning the partition function (see Eqs. flA.ip and (IA.2p i such that 
the exponential of the source terms is replaced by the constant and linear terms in the 
power series expansion 

A„( A, /i, A V) = n trs. *■' 


y 

X J](1 + /iJi[5,] + /3J3[^.]). (A.16) 
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The structure in Eq. lA. 161 ensures that there will be at most one single source term per 
lattice site and otherwise does not have any impact on the derivation presented above. 
From here on, we will only consider the partition function presented in Eq. lA. 161 for the 
branched polymers system. 


Mean Field Hamiltonian 


We can now convert the lattice model over a discrete set of S'- vectors into a continuous 
held theory 'if{x) using a Hubbard-Stratonovich transformation and connect the lattice 
fugacities K, fi, and /a to physical quantities of chemical potential and concentration 
in the mean held approximation. 

It is necessary to carefully treat the sum over nearest neighbors on the lattice given 
in Eq. flA.16|) in order to change the lattice model into a continuous held theory. The 
sum can be written as a double sum over lattice sites multiplied by a nearest neighbor 
delta function 


{x,x') 


V V 

9 • 


(A.17) 


The function 5(^x,x') is similar to a Kronecker delta function, and is explicitly dehned 
as an operator that evaluates to 1 when x and x' are neighbors and 0 otherwise. The 
additional factor of half prevents double counting. Using Eq. flA.17j) . we now perform a 
Hubbard-Stratonovich transformation to introduce the auxiliary held 'ip{x) 


■y E 

g x,x' 


'■oo E s.f ,.-tp{xy'4>{x')+j2VKit’ixySx 

Vfj e " , (A.18) 


where is the inverse of the Si^x,x') operator. Using Eq. (lA.lSp . the partition function 
Eq. flA.16p can be written as 


Z^{K,h,h-V)= \ Vfje 






X 


l[tTs^ + /i Ji[5,] + fMS,]). (A.19) 


The hrst term in the second line of Eq. flA.lhp is simply the generating function 
performed in Eq. flA.9p . We now use Eqs. flA.ldaP and flA.13bp to evaluate the integrals 
associated with the sources Ji and J3 in Eq. flA.19p in the n —)■ 0 limit. Without loss of 
generality, the source terms can pick a special direction. To simplify the integrations, 
we choose the (1, 0,..., 0) direction and thus Eqs. flA.13aP and flA.13bp can be written 


limtrs = dn, 


n —>0 


hmtrsS'iS'jS'fcJ3[5] = dudjidki- 
rn-O 


(A.20) 

(A.21) 
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Inserting Eqs. flA.20p and flA.2ip into Eq. (IA.19P and performing the integral over S- 
vector, the partition fnnction in the n —)■ 0 limit becomes 


Zn{KJ,J,-V)= /P^e 




V 


^ n ( ^ + ylV^(a^)r + /iV^V'i(a;) + ). (A.22) 


The fifji^x) and f^^jf^x) terms are proportional to the end and branch-point densities, 
while '0^ is proportional to the monomer density. Since for most physically relevant 
systems the ratio of branch or end-points to monomers is low, the fi'ijji{x) and fsi^fix) 
terms will be mnch smaller than the 'ijj‘^{x) term. Raising the second line of Eq. flA.22p 
into the exponent (1 -|- X = and expanding the logarithm, we dehne a new 

effective Hamiltonian as given in Eq. (jSp. 
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